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3.5 Karman's momentum integral approach 

Ref: Schlichting: Boundary layer theory 

With a general pressure gradient the boundary layer equations can be solved by a va- 
riety of modern numerical means. An alternative which can still be employed to simplify 
calculations is the momentum integral method of Karman. We explain this method for a 
transient boundary layer along the x-axis forced by an unsteady pressure gradient outside. 
This pressure gradient can be due to some unsteady and nonuniform flow such as waves or 
gust. 

3.5.1 General ideas 

Let us limit our discussoin to a two dimentional flow. Recall the continuity 

u x + w z = (3.5.1) 

and boundary layer equation: 

p(u t + uu x + wu z ) = —p x + pu zz (3.5.2) 
= - Pz (3.5.3) 

The boundary conditions are 

u — w — 0, z = 0; u — > U, w — > 0, z — > oo. (3.5.4) 
For general U (x, t) the longitudinal pressure gradient is 

-p x = p(U t + UU X ) (3.5.5) 

so that 

p(u t + uu x + wu z ) = p(U t + UU X ) + pu zz (3.5.6) 

Instead of solving the initial-boundary-value problem accurately we introduce a moment 
method by integrating the momentum equation across the boundary layer, then make a rea- 
sonable assumption on the velocity profile to get the longitudinal variation of the boundary 
layer thickness. This is called the Karman's momentum integral approximation in boundary 
layer theory. 



Let us add the following two equations 

—fiu zz = p{U t + UU X - (ut + uu x + wu z )} 



and 

to get, 



= (U — u)u x + (U — u)w z 



-»u zz = p fJL(U-u) + ^[u(U - «)] + (U- u) 9 ^ + £-[(w(U - «)]} 

Now let us integrate across the boundary layer, 

{d f°° d f°° dU f°° 1 

— y o (u - u )dz + — y o u (u - u)dz + — y o (u -u)dz\ + P [w(u -u)\™ 

The boundary terms again vanish. Finally, we have, 



-J o p(U-u)d Z + -j Q P (Uu-u)dz+—j o (U-u)dz=p- 



d 



dU 



du 



(3.5.7) 



This is the Karman momentum integral equation, representing the momentum balance across 
the thickness of the boundary layer. 




Figure 3.5.1: Displacement thickness 
Karman also introduces the displacement thickness as a measure of the lost volume 

[CO [CO / fl,\ 

U5 1 = J (U- u) dz, or <?! = J ( 1 — - J cfe, (3.5.8) 
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and the momentum thickness as a measure of the lost momentum 



u 2 5> 



Jo 



/ u(U — u)dz, or 5- 



(3.5.9) 



both due to the slowing down of the fluid in the boundary layer. See Figure 3.5.1. 

In certain cases when the velocity profile in the boundary layer can be reasonably guessed 
in advance, the Karman momentum intergal equation can be the basis of obtaining an 
approximate solution. The procedure is to assume a reasonable profile : 



and substitute (3.5.10 ) into (3.5.7) equation. After evaluating the integrals a differential 
equation is obtained for the boundary layer thickness 5(x, t), which can be more easily solved 
for certain initial and boundary conditions. 

3.5.2 Application to a transient boundary layer along a flat plate 

Consider the two dimensional boundary layer near edge of half infinite plane along the x 
axis due to the impulsive start along its own plane. There is no motion anywhere for t < 0. 
At t — the plane suddenly advances from right to left perpendicular to its leading edge. 
What is the boundary layer flow? Referring to Figure 3.5.2 where the coordinate system is 




(3.5.10) 



U=0 




Figure 3.5.2: Front of boundary layer 



fixed on the plane. For allt > 0, 




dU _ dU 
dt dx 



The boundary layer equation reads: 



p{u t + uu x + wu z ) = jJLU 



'ZZ 



(3.5.11) 



and the Karman momentum inetgral equation reads 




(3.5.12) 
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Let us assume 



then 



and 



/•oo 

5 1 = 5 (1 - /(?7)W?7 = ai<J (3.5.14) 
J o 

/•oo 

<5 2 = 5/ f(r,)(l-f(r,))dr, = a 2 8 (3.5.15) 

u, = ^/'(0). (3.5.16) 
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Substituting into Eq. (3.5.12) 
Th.cr6for6 

Mffl + ^m^^O) ,3.5,7) 

which is a first order wave (hyperbolic) equation for 5 2 /2. We must add the boundary and 
initial conditions: 

<?(0,t) = 0, Vt>0 (3.5.18) 

5(x,0) = 0, \/x>0 (3.5.19) 

Heuristic solution 

The fact that the boundary value in (3.5.18) is independent of t for all t suggests that the 
solution is independent t for sufficiently large x, i.e., 

— — = I \ ; 3.5.20 

dx \ 2 J a 2 U v ; 

subject to the boundary condition (3.5.18). The solution is simply 



This is the approximate version of the solution by Blasius who solved the steady boundary 
layer equation 

uu x + vu y = vuy y Ioy x > , 1) > (3.5.22) 

by the method of similarity. 

On the other hand, the fact that the initial value in (3.5.18) is independent of x for all x 
suggests that the solution is independent of x for sufficiently large t, i.e., 
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subject to the initial condition (3.5.19). The solution is simply 



S =,P» (3.5.24) 

This is the approximate version of the solution by Rayleigh's problem which is known to be 
governed by 

u t = uuyy for t > 0, y > (3.5.25) 
To find the range of each solution in the x ~ t plane, we equate the two <5's and get 

x = ^ (3.5.26) 

Thus (3.5.21) holds in the wedge x > ^ > and (3.5.24) holds in the wedge < x < 
See Figure 3.5.3. 

Physically: if t > aix/a 2 U, (3.5.21) applies and one has the steady Blasius flow past a 
semi-infinite plate; the boundary layer is already in the steady state. See Figure 3.5.4, If 
t < aix/a 2 U, (3.5.24) applies and one has Rayleigh's problem of impulsively started plane 
infinite in both x < and x > 0. The boundary layer is still in the initial stage and the 
effect of the leading edge is not felt elsewhere. The result is summarized in Figure 3.5.5. 

To calculate ot\ and a 2 let us make a special choice of the velocity profile so that /(0) = 
0, /'(I) = /"(I) = (due to Karman and Polhausen) 

f( v ) = 2r] - 2rf + r] A , 0<r]<l (3.5.27) 



Note that 
then 



/(0) =0, f{r j ) = 2-6r7 2 + V, /" (rj) = -12r) + 12 V (3.5.28) 

ai = ^ (l - 2r] + 2if - rf) dr] = 3/10 = 0.3 

a 2 = J (277 — 2?7 3 + ?7 4 ) (l-2?7 + 2r] 3 - r] 4 ) dr] = 37/315 = 0.117 



and 



m = % or I 



2U 



z=0 

Hence the displacement thickness of the Blasius boundary layer is 



6 1 = ai 6 = a J*^ = i™M (3.5.29) 

V 012U v u 

The momentum integral method is the special case of the moment method, since the 
Karman equation is the zeroth moment of the boundary layer equation. For the classical 
steady boundary layer problem solved exactly by Blasius using the similarity method, the 
momentum integral approximation gives fairly good results, even with various crude profiles, 
see Table 3.1. 
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Figure 3.5.3: Evolution of a transient boundary layer, (a) Blasius region and Rayleigh region 

.5 -t 




Figure 3.5.4: Evolution of a transient boundary layer, (b) Blasius region 
Formal solution 

A more formal solution can be obtained as follows: Let 

a — a.\ b = a 2 U c = vf(0) 

then 

d f5 2 \ , d f5 2 \ 



(3.5.30) 



The solution can be facilitated by a change of coordinates from (x, t) to (£, () where 



then 



£ = ax + bt ( = ax — bt 

d__d_d£ d_d(_ d_ d_d^_ d_d^ 

dt d£ dt d( dt ' dx d£ dx d( dx 

(6% = (««),«.+ (*•),« 



(3.5.31) 



Therefore, from (3.5.17), 

1 

2 



ab(5 2 ) -ab(f) +ab(5 2 ) +ab(& 
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Figure 3.5.5: Evolution of a transient boundary layer, (c) Rayleigh region 



Velocty profile 


Displacment boundary layer thickness 




W- 

1 V vx 


7] 


1.732 


sinf 


1.741 


2r] - 2 rf + r] 4 


1.754 


Blausius 


1.721 



Table 3.1: Comparson of approximate solution by momentum integral method with the exact 
solution of Blasius 



Integrate once 



2 V" Jt 2ab 2aia 2 U 



(3.5.32) 



2 2a x a 2 U 
where G is an arbitrary function of (. 

To determine G(() we first use the initial condition that 5 = at t = for all x > . 

0=i£*°L(£)„ + G(O« 



or 



2aia2^ 



2a 1 a 2 £7 



ax) + G(ax) for x > 0. 



Therefore 



B(0 = -^, <>° 



(3.5.33) 



2aiQ;2?7 

What is G(C) for ( < 0? Let us use the boundary condition at x = that 5 = for all 
t > 0. 

*/'(0) 



or 



= 



2a.\OL 2 U 



(0«=o + G(C)x=o 



2a\a 2 U 



(-bt) +G(-bt), t>0 
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Note that for t > 0, -bt < 0. Therefore, 

G(C) = ^^ 7 C for all C<0. (3.5.34) 

Eqs. (3.5.33) and (3.5.34) complete the solution for all ( > and ( < 0. 
Let us return to x and t 

— = —£ - —Q = —2bt for ax - bt > 

2 2a\a 2 u 2a\a2<J 2aia2U 

4 + - — C — n 7720a; for ax — bt <0 



2aia 2 U 2aia 2 U 2a\a 2 U 

Therefore, 



' 2 "/'(°>. ( for x> ^ t (3.5.35) 



and 



*=,H^z for x<^ ( (3.5.36) 

V Oi 2 U Cti 



These are the same as (3.5.24) and (3.5.21). 



